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ABSTRACT 

We present results of a gravitational-lensing and optical study of MACS J1423. 8+2404 
(z=0.545, MACSJ1423), the most relaxed cluster in the high-redshift subsample of 
clusters discovered in the MAssive Cluster Survey (MACS). Our analysis uses high- 
resolution images taken with the Hubble Space Telescope in the F555W and F814W 
passbands, ground based imaging in eight optical and near-infrared filters obtained 
with Subaru and CFHT, as well as extensive spectroscopic data gathered with the 
Keck telescopes. At optical wavelengths the cluster exhibits no sign of substructure 
and is dominated by a cD galaxy that is 2.1 magnitudes (K-band) brighter than the 
second brightest cluster member, suggesting that MACS J1423 is close to be fully viri- 
alized. Analysis of the redshift distribution of 140 cluster members reveals a Gaussian 
distribution, mildly disturbed by the presence of a loose galaxy group that may be 
falling into the cluster along the line of sight. Combining strong- lensing constraints 
from two spectroscopically confirmed multiple-image systems near the cluster core 
with a weak- lensing measurement of the gravitational shear on larger scales, we derive 
a parametric mass model for the mass distribution. All constraints can be satisfied 
by a uni-modal mass distribution centred on the cD galaxy and exhibiting very little 
substructure. The derived projected mass of M(< 65"[415kpc]) = (4.3 ± 0.6) x 10 14 
Mq is about 30% higher than the one derived from X-ray analyses assuming spherical 
symmetry, suggesting a slightly prolate mass distribution consistent with the opti- 
cal indication of residual line-of-sight structure. The similarity in shape and excellent 
alignment of the centroids of the total mass, K-band light, and intra-cluster gas dis- 
tributions add to the picture of a highly evolved system. The existence of a massive 
cluster like MACS J1423, nearly fully virialized only ~ 7 Gyr after the Big Bang, 
may have important implications for models of structure formation and evolution on 
cosmological time scales. 

Key words: Gravitational lensing: strong lensing - Galaxies: clusters: individual 
(MACS J1423.8+2404) 



1 INTRODUCTION 



The MAssive Cluster Survey (MACS, [Ebeling et al.pOOl ) 
has compiled a complete sample of 12 very X-ray luminous 
galaxy clusters at z>0.5 ( |Ebeling et al.|[2007 ), providing a 
unique opportunity for comprehensive studies of the densest 



regions of the cosmic web at intermediate redshifts. As noted 



by Kartaltepe et al. (2008), MACS J1423 stands out among 



these 12 extreme systems as the most dynamically relaxed; 
it is, indeed, the most massive cool-core cluster known at 
these redshifts. The system has been used in several sta- 
tistical studies ( LaRoque et al.|[2003 Puetzfeld et al.| 2005 
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Figure 1. Left: MACS J1423 as observed with HST/ACS (152" xl52", corresponding to 965x965 kpc 2 ). Multiple-image systems are 
labelled, the weak-shear field is overlaid in cyan. Right: The top four panels (4.7 /; x4.7 /; , 30x30 kpc 2 ) show system 1-2, assumed to 
belong to the same background source galaxy. System 1 is systematically brighter than system 2. The fifth panel (7. 9" x 7.9", 50x50 
kpc 2 ) shows image 1.5 after the light of the cD galaxy has been subtracted. The cyan circle highlights regions of star formation commonly 
found in the cD galaxies of cool-core clusters. Indeed, the cD shows strong line emission in IFU observations (Edge et al., in prep.), 
indicative of ongoing cooling. The final three panels show images 3.1, 3.2 and 3.3 respectively (9.4" x9.4", 60x60 kpc 2 ). 



Bonamente et al.|[2006| |Kotov &; Vik hlinin 2006} where its 
relaxed dynamical state has been particularly important for 



cosmological work using the baryon fraction ( 


Allen et al. 


2004| ISchmidt & Allen||2007| |Allen et al.||2008| 


Ettori et al. 


2009). Since the abundance of massive relaxed clusters at 



high redshift has important implications for theoretical and 
numerical models of structure formation and evolution, we 
here attempt to further test the system's relaxation state 
using gravitational lensing. All our results use the ACDM 
concordance cosmology with = 0.3, Ha = 0.7, and a 
Hubble constant Ho = 70 kins" 1 Mpc -1 . Magnitudes are 
quoted in the AB system. 



2 OBSERVATIONS 
2.1 Imaging 

Ground based panoramic imaging of MACS J 1423 was per- 
formed in the B, V, R c , I c , and z bands with the Suprime- 
Cam camera on the Subaru telescope, and in the u* and 
K bands with MEGACAM and WIRCAM on CFHT. The 
resulting object catalogues were used to compute photomet- 
ric redshifts for all galaxies in a 0.5 x 0.5 deg 2 field follow- 
ing the methodology described in |Ma et al.| (|2008). Cluster 



members were defined to be galaxies with photometric or 
spectroscopic redshifts (see below) within ±0.05 of the clus- 
ter redshift. 

MACS J1423 was observed on 2004, June 16 with the 
Advanced Camera for Surveys (ACS) on-board HST for 4.5 
ks and 4.6 ks through the F555W and F814W filters, respec- 
tively (GO: 9722, PI: Ebeling). The data were reduced using 
an automated pipeline developed by the HAGGLeS team 
(Marshall et al., in prep.), that uses multidrizzle version 
2.7.0. Following manual masking of satellite trails, cosmic 
ray clusters, and scattered light in the flat-fielded frames, 
all exposures were carefully registered onto a common 0.03 
pixel grid ( Schrabback et al.|2007| , generating the shift files 
needed by multidrizzle. The individual frames in each fil- 
ter were then drizzled together using a square kernel with 
pixfrac 0.8, applying updated bad pixel masks and optimal 
weights ( |Schrabback et aL][2 009). 

In spite of the lensing efficiency decreasing with in- 
creasing cluster redshift, MACS J1423 exhibits clear strong- 
lensing features. Visual inspection of the ACS frames imme- 
diately reveals two multiple-image systems (Fig.[T]). The ob- 
served configuration of lensing features for the first of these 
can be explained by a single background galaxy, resolved 
into two components, each of which represents a multiple- 
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Figure 2. Co-added spectra for system 1-2 at z = 2.84 (top) and 
system 3 at z = 1.779 (bottom). 



image system (referred to as 1-2 in the following) in a typi- 
cal Einstein-cross configuration. The mass model predicts a 
central fifth image which is detected at the predicted posi- 
tion under the cD galaxy whose light has been subtracted 
using the IRAF routine ellipse. The other multiply imaged 
system, #3, constitutes a typical naked-cusp configuration, 
with three images on the same side of the cluster. 



2.2 Arc spectroscopy 

Spectroscopic observations of all lensing features were per- 
formed on March 6, 2005 using the Low Resolution Imager 
and Spectrograph (LRIS, |Oke et al.||1995) of the Keck I 
telescope in multi-object spectroscopy mode. We used the 
400/3400 grism, D560 dichroic, and 400/8500 grating cen- 
tred at 6800A to achieve continuous wavelength coverage 
from 3200 to 8000A. The reduction of the data obtained in a 
total integration time of 10.8 ks followed the standard steps 
of bias subtraction, slit identification and extraction, flat 
fielding, wavelength calibration, and flux calibration, and 
was performed using the python routines described by |Kel-| 
son et al. (2002). From the spectra thus obtained we mea- 



sured redshifts of z=2.84 for the images labelled 1-2.1, 1-2.2 
and 1-2.3; the redshift of features 3.1 and 3.2 was found to be 
z = 1.779. The spectra for each component show identical 
features (Lymana in emission for system 1-2, UV absorp- 
tion lines in system 3), therefore they have been co-added 
on Fig. [2] 



2.3 Galaxy spectroscopy 

Using the LRIS and DEIMOS spectrographs on the twin 
Keck- 10m telescopes on Mauna Kea we conducted an ex- 
tensive spectroscopic survey of galaxies in the field of 
MACS J 1423. Redshifts were measured for 396 objects, 60 of 



them within the solid angle covered by our HST/ACS obser- 
vation. All redshifts were used to calibrate the photometric 
redshifts mentioned above. Using iterative 3a clipping we 
find z c i=0.5450 and cr=1320lg^ km s _1 for the cluster red- 
shift and velocity dispersion, respectively (n z =140, Fig. |3j). 
When limited to the ACS field of view, the same analysis 



yields z c i=0.5424 and cr=1470] 



km 



=40). Even 



for a system as X -ray luminous as M ACSJ1423 (Lx,boi = 
3.7 x 10 45 erg s~ 1 , |Ebeling et al.|2007) , t his velocity disper- 



sion falls high on the Lx-cr relation (e.g. Mahdavi & Geller 



2001), without, however, constituting a significant outlier. 

Inspection of the galaxy velocity histogram (Fig. [3] top) 
shows the distribution to be mildly skewed toward lower red- 
shifts; a one-sided Kolmogorov-Smirnov test finds it to be 
inconsistent with a Gaussian at the 2a confidence level. We 
investigate whether the two peaks at z^0.531 and z^0.553 
represent subclusters along the line of sight by selecting 
three galaxy subsamples as indicated by the coloured bands. 
The distribution of the respective galaxies on the sky is 
shown in Fig. [3] (bottom). Unlike the galaxies with 2^0.553, 
we find the galaxies in the peak around z^0.531 to be sig- 
nificantly (3.3<j) more concentrated in projection on the 
sky than the overall cluster galaxy population, suggesting 
a gravitationally bound substructure viewed almost exactly 
along the line of sight. We thus conclude that the velocity 
dispersion of MACS J1423 is mildly inflated by a possible 
infall of a loose group at z = 0.531, viewed in projection 
and superimposed directly onto the cluster core. 



3 DISTRIBUTION OF LIGHT & MASS 
3.1 Light Distribution 

We use the K-band magnitudes of (photometric) cluster 
members to create a luminosity-density map of the cluster 
light. By fitting a 2D Sersic profile to the cD galaxy which 
dominates the K-band cluster light (we use the GALFIT soft- 



ware of |Peng et al. (2002)), we measure an integrated mag- 



nitude of K=16.27 for the cD, compared to K=18.33 for the 
second brightest cluster member. The magnitude difference 



is K g 



=2.06. 



3.2 Lensing Analysis 

In addition to the two spectroscopically confirmed strong- 
lensing features discussed above, we also take into account 
the weak- lensing signal as reflected in the sheared F814W 
images of the background galaxy population. Details of 
the weak-lensing analysis will be presented in a forthcom- 
ing publication (Jauzac et al., in prep.). In short, we fol- 
low the weak lensing pipeline developped by Leauthaud 
et al. ( 2007| 2010), but adapted the code to work on clus- 
ter field, in order to efficiently remove the galaxy cluster 
members and the foreground galaxies that would other- 
wise reduce the strength of the shear signal. In doing so, 
we consider only background sources located outside the 
Einstein radius of, for MACS J1423, R E = 21" (here we 
refer to the effective Einstein for a source at z = 2, de- 
fined as the angular radius from the centre of the cluster 
at which the average convergence equals 1), thereby min- 
imising any systematic bias from the strong-shear regime 
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Figure 3. Top: Redshift histogram of galaxies in the field of 
MACSJ1423. The dashed line shows the best-fitting Gaussian 
(z cl =0.5450; cr=1320tgi km s _1 , n z =U0). The vertical dotted 
lines mark the zblcr range; the coloured bands are designed to sep- 
arate radial velocities of potential fore- or background structure. 
Bottom: The distribution of selected cluster members on the sky 
(same colour scheme as in top figure: the main cluster galaxies 
are shown as green circles, and the populations of the two peaks 
are shown as filled diamonds). Galaxies with z ~ 0.531 (blue) are 
tightly concentrated and thus likely to constitute a bound galaxy 
group, falling into the cluster core along our line of sight. The 
black square outlines the area covered by our HST/ACS observa- 
tion. The dashed lines delineate small areas that are not covered 
by our spectroscopic survey. The small black dots represent other 
galaxies with spectroscopic redshifts but outside the cluster. 



(Mas sey fe Goldbe rg 2008). The 192 galaxies thus selected 
within the ACS field (corresponding to a number density of 
16 galaxies per square arcminute) represent 192 constraints 
in the weak-lensing analysis. Adding the 18 constraints from 
the two strongly lensed multiple-image systems yields a to- 
tal of 210 constraints, all of which are used in order to 
build a parametric mass distribution within the framework 



of the publicly available lenstool software (| Julio et al. 
2007 http://www.oamp.fr/cosmology/lenstool/) in which 
we have implemented the inclusion of the weak lensing con- 
straints. The optimisation is performed in the image plane. 

To characterise the mass distribution in MACS J 1423 
we adopt a dual Pseudo Isothermal Elliptical Mass Distri- 



bution (dPIE, |Limousin et al.|2005| |EKasd6ttir et al.|2007| ), 
parametrised by a fiducial velocity dispersion cr, a core ra- 
dius r CO re and a scale radius r s (set to 1000"). In addition 
to this large-scale mass component, which models both the 
intracluster gas and the far larger dark matter component 



( Bra dac et al.j 2008 ; Sommer-Larsen & Limousin 2009 ) , we 
also take into account the perturbations associated with in- 
dividual cluster galaxies. Specifically, we include in the mod- 
elling all cluster members brighter than 22.5 (F814) and 
within a projected cluster-centric distance of 35" from the 
cD galaxy. Empirical scaling relations (without any scatter) 
are used to relate their dynamical dPIE parameters (central 
velocity dispersion and scale radius) to their luminosity (the 
core radius being set to 0), whereas all geometrical param- 
eters (centre, ellipticity, position angle) are set to the val- 
ues measured from the light distribution (see, e.g. Limousin 
et al.|2007b Richard et al.|2009a| . In total our mass model 



is thus described by eight free parameters. The range of al- 
lowed parameter values characterising the dPIE mass com- 
ponent is limited by priors as follows. Position: ± 25" along 
the X and Y directions; ellipticity: <0.6 (expressed in units 
(a 2 — b 2 )/(a 2 + 6 2 ), where a and b refers to the semi-major 
and semi-minor axes of the ellipse describing the mass distri- 
bution); position angle: 0-180 degrees; velocity dispersion: 
500-1500 krns" 1 ; core radius: 1-20". For the galaxy-scale 
component, we allow the velocity dispersion to vary between 
100 and 250 krns -1 , whereas the scale radius was forced to 
be less than 70 kpc in order to account for tidal stripping 
of their dark matter haloes (see, e.g. |Limousin et al.|2007a 



2009 1 |Natarajan et al.|2009| fWetzel fe White||2009[ and ref- 



erences therein. 



3.3 A Uni-modal Mass Distribution 

The shear field constructed from our analysis of the images 
of the weakly lensed background galaxy population (Fig. [I]) 
already suggests a uni-modal mass distribution centred on 
the cD galaxy. Combining strong- and weak-lensing con- 
straints, we confirm that a single cluster-scale mass com- 
ponent is able to reproduce the full set of observational 
constraints, with a reduced x 2 value of 175/202. The cor- 
responding best-fit values of all model parameters are sum- 
marised in Table 1. We find the centre of the elliptical dPIE 
mass distribution to coincide with that of the cD galaxy; its 
position angle of 116=b2 degrees matches the orientation of 
the cD galaxy (126=bl0 degrees). No additional large-scale 
mass component is required by the data. 

Fig. [4] compares the 2D mass map corresponding to the 
best-fit lens model to the X-ray surface brightness as ob- 
served with Chandra (from Ebeling et al. (2007)) and to the 
K-band light map, both adaptively smoothed to 3a signifi- 
cance using ASMOOTH (Eb eling et al.|2006| . The agreement 
between these independent observational tracers of the three 
cluster constituents (dark matter, gas, and galaxies) is re- 
markable. 



4 DISCUSSION & CONCLUSIONS 

Using optical and lensing observations of MACS J 1423 we 
find overwhelming evidence for a uni-modal mass distribu- 
tion whose centre, ellipticity, and orientation are in excel- 
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Component 5x{") 5y(") e Or (kpc) a (kms 1 ) 

Cluster Halo -1.0±1.2 O.Oil.O 0.4±0.1 116.0±1.7 78.7±25.2 1081±110 
Perturbers - 32±|£ 14ljl|f 

Table 1. Best-fit values of all mass-model parameters. Coordinates are defined with respect to the cD galaxy. Error bars correspond to 
3cr confidence level as inferred from the MCMC optimisation. For the dPIE mass component, r corresponds to the core radius whereas for 
the galaxy-scale perturbations it corresponds to the scale radius. 




fidence level, Schmidt & Allen 2007} from an analysis of 



Right Ascension (2000) 



Figure 4. The distribution of mass, gas, and galaxies in the core 
of MACS J 1423, as reflected by the adaptively smoothed K-band 
luminosity of cluster galaxies (red), the total mass derived from 
the gravitational-lens model (black) and the adaptively smoothed 
X-ray surface brightness from Chandra observations (blue). 



lent agreement with the corresponding parameters of the 
cluster light and intra-cluster gas distributions (Fig. [4]). We 
quantify the dominance of a single, central mass compo- 
nent by computing the substructure fraction, / su b, defined 
as the mass within a given cluster-centric radius that can 
be attributed to substructure, divided by the total mass 
within the same radius (ISmith & Taylor||2008| . In the case 
of MACS J1423, only cluster galaxies (the cD galaxy is by 
definition excluded) contribute to the substructure fraction. 
We find / su b(< 250 kpc)=0.048±0.025, comparable to values 
measured for the most relaxed clusters at z ~ 0.2 ( |Richard| 
et al.||2009b| . Within the high-redshift MACS subsample, 
the substructure fraction of f S ub(< 500 kpc)=0.045=L0.025 
of MACS J1423 is in stark contrast to the value of 0.25±0.12 
measured for MACS J1149. 5+2223, the most complex clus- 
ter lens studied to date (Sm ith et al.| 2009). 

The radial mass profile of MACS J1423 as derived from 
our lensing analysis is well constrained across the full ACS 
field of view, i.e., out to a radius of 100" (640 kpc). Previous 



result s for the mass within 65" (415 kpc) of 5.0+o g x 10 



the Chandra X-ray data, agree within the errors with our 
measurement of 4.35 ± 0.6 x 10 14 M (3a). We note, how- 
ever, that the existing X-ray analyses of MACS J1423 as- 
sume spherical symmetry which will bias the resulting mass 
estimates low for a prolate matter distribution. The differ- 
ence of rsj 30% between the lensing and X-ray mass may 
thus be due to an elongation of the system along the line of 
sight, consistent with our finding of residual accretion onto 
the cluster core (Section 2.3). 

Richard et al. (2009b) investigated the relationship be- 



tween the effective Einstein radius and the strong lensing 
mass (measured within 250 kpc) for a sample of 20 strong 
lensing clusters at z ~ 0.2 drawn from the Local Cluster 
Substructure Survey (LoCuSS). According to their statisti- 
cal analysis, one would expect a mass of 2.55^; 9 6 ° 5 10 14 M 





for an undisturbed cluster with an Einstein radius of 21". 
The excellent agreement with our measurement of 2. 46 ±0.2 
10 14 Mq underlines again the relaxed dynamical state of 
MACS J1423. 

Although our analysis of the radial velocity distribu- 
tion of cluster members suggests that MACS J1423 evolved 



through mergers along a line-of-sight filament (Section 2.3), 
the rate of current or recent infall is low and does not result 
in a noticeable perturbation of the cluster core. In fact, the 
large luminosity gap of K gap =2.06 mag between the cD and 
the second brightest cluster galaxy, as well as the excellent 
match of the dark-matter, intra-cluster gas, and cluster light 
distributions, represent strong evidence for a dynamically 
highly evolved cluster core ( | Jones et al.|2 003 Milosavljevic 



et al. 2006), undisturbed by major mergers for at least 1 



M , ( |LaRoque et al.|20 03P| using observations of the Sun- 
yaev-Zel'dovich effect, and of 3.1±S;? x 10 14 M (3a con- 



Gyr, the characteristic time scale for relaxation via dynam- 
ical friction. The existence of MACS J1423 only 7 Gyr after 
the Big Bang thus adds an important observational con- 
straint to N-body simulations of structure formation and 
evolution in the current cosmological model. 
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